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PREFACE 


This investigation of stability analysis of the atmosphere was 
conducted at the U. S. Naval Postgraduate School, Annapolis, dryland 
during the -winter and spring of 1948 in partial fulfillment of the 
requirements for the degree of Master of Science in Aerology. The 
object of the work was to check the findings of previous investigations 
along this line by taking observations over a longer period of time and 
to determine the desirability of using the Slice Method as a means of 
objective forecasting using single station soundings. 

The writer is indebted to the staff of the U. S. Naval Post¬ 
graduate School for their guidance, and especially to Professor W. D. 
Duthie whose helpful suggestions made possible a more thorough investi¬ 
gation into the field of stability analysis. The cheerful response of 
the U. S. Weather Bureau to requests for information and data made the 
initial steps in the project much simpler. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


C 

Po 

p, 

To 

T, 

<T 

Xfn 

R 

3 

T* 

cr'&t 


Jfe') 



A, 


Absolute circulation acceleration. 

Pressure at the bottom of a finite layer of the atmosphere. 
Pressure at the top of a finite layer of the atmosphere. 
Temperature at the bottom of a finite layer of the atmosphere. 
Temperature at the top of a finite layer of the atmosphere. 
Thickness of a finite layer of the atmosphere. 

Actual lapse rate of temperature -with height. 

Dry adiabatic lapse rate. 

Moist adiabatic lapse rate. 

The gas constant for air. 

The acceleration of gravity. 

Mean virtual temperature of a layer of the atmosphere. 

A vertical impulse, or velocity pertubation exerted over a 
time & t. Assumed arbitrarily to be unity in the Slice Method. 
The ratio of the mass of ascending air to the mass of descending 
air in a convective cell. 

A dimensionless numerical function giving a measure, of 
stability of a convective cell with dry ascent dry descent. 

A dimensionless numerical function giving a measure of the 
stability of a convective cell having moist ascent-dry descent. 

A dimensionless numerical function giving a measure of stability 
of a connective cell having moist ascent- moist descent. 

A dimensionless numerical function, connoting difference, which 
when added to D^) gives W (tr) . 

Same as above except that when added to D(*s) gives II(^) . 
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TABLE. OF SYMBOLS AND ABBREVIATIONS 
(continued) 


mbs. Millibars, 

MCCL Mean convective condensation level. 



CHAPTER I 


INTRODUCTION 

Stability analysis by the Slice, or Circulation acceleration, 
method has received little use outside of the U. S. Naval Postgraduate 
School. The original development by Beers of the msthod was left in a 
rather cumbersome form for computing stability, so that the work 
invo3.ved in making computations in many cases ms prohibitive. 

Duthie, with no further assumptions, reduced the basic equations 
to a very simple computing form. 

Next came the problem of testing the method in a statistical study. 
Preliminary investigation along this line was conducted by Cain, Ashley, 
and Champion, in a seminar at the Postgraduate School. Their .results 
were good, and indicated that there ms merit in this method of stability 
analysis. 

Cain, Ashley, and Champion devoted most of their work to investi¬ 
gating stability over large areas of the weather map. They found, in 
general, that isolines of stability (or instability) could be drawn such 
that most of the thunderstorms and convective activity occurred within 
regions of instability, i.e., positive values of Circulation Acceleration 
(hereafter designated C for brevity) and that these regions of instability 
could be followed from map to map. They also investigated the correlation 
of convective activity with instability for a single station (Washington, 
D. C.) for a period of several days. 

The present work was an exbention of the latter. All atmospheric 
soundings for Washington, D. C. for the months of July, and August 1947 
and the associated weather were obtained from original data from the 


1 



. M' 

U. S. Weather Bureau. The values of C were computed for || of .2 and 

.3, from the Mean Convective Condensation Level (MCCL) to 500 mb., 600 

mb. and the 0° isotherm. Various synoptic correlations were investi- 

gated, i.e., C, moisture (mixing ratio, and Relative humidity), mean 

o 

temperature of the 1000 mb.-700 mb. layer, height of 0 isotherm, 
temperature at the 700 mb. surface, and distribution of moisture. C 
was computed for the lifted and non-lifted soundings prior to a squall- 
line passage or frontal passage. Since some unstable values of C did 

not give thunderstorms and some thunderstorms occurred with stable 

* 

values of C, various weight factors were applied in an attempt to make 
all data consistent. Moist ascent - moist descent was applied under 
certain conditions._ 

It should be noted that except where indicated the correlation is 

very nearly synoptic. Both 1100 and 2300 EST soundings, and U. S. 

Weather Bureau analysed maps of 1330 and 0130 EST were used. 

This investigation although not producing spectacular results, 

indicates the lines along which further research should be conducted, 

and points out some fallacies in previous methods of computation. The 
o . 

use of 70 /o mean relative humidity as the dividing line between wet 

and dry layers is definitely not desirable. Rather, the degree of 

saturation of the atmosphere in the lower layers, specifically, that 

below 700 millibars seems to be more important in determining the prob- 

ability of occurrence of thunderstorms. The computation of C to the 600 

millibar level appears to produce a better measure of stability and of 

the liklihood of occurrence of thunderstorms than the 0° C. level or the 

500 millibar level. The value of jf of seems to give more accuracy 
3 

than •jjjj when using the number zero as the reference level for positive 
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and negative stability, however this is open to question. The 2300 
sounding proved to have very little value for forecasting purposes. The 
1100 sounding, on the other hand, gave very good results for short range 
forecasts* The practice of lifting the sounding 100 millibars for an 
expected cold frontal passage appears sound and by doing this the 
forecasting accuracy was improved considerably. 
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CHAPTER II 


EXPLANATION OF THE SLICE METHOD 

The Slice Method of stability analysis as originally presented 
by Bjerknes (1) and adapted by Petterssen (2) classifies the stability 
of a single level of the atmosphere of unit thickness. Beers (3) 
modified the method so that account may be taken of the stability layer 
by layer through the atmosphere, especially through the column where 
cumulus activity may occur, the layers being of arbitrary thickness. 

His development is not presented here, but his results and their modifi¬ 
cations by Duthie (4) are interesting and mil be shown. 

Beers arrived at three essentially different equations respecting 
stability, depending upon the moisture distribution in the atmosphere. 
These cases and their appropriate equations follow. When the Slice 
Method is referred to hereafter Beers' treatment will be implied. 

I. Dry Layers. Both ascending and descending currents are 
initially and remain so dry that condensation does not occur. Then 

(l) c - R /o [(*- *)(' + 

II. Liquid Water Present. The layer contains enough liquid 
water so that both ascending and descending currents follow the saturated 
adiabat, i.e., condensation occurs in the ascending current while evapo¬ 
ration occurs in the descending current. Then 

(?.) c - R In [(ir- r^)(i + *'&t. 
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III. Wet Layers. The ascending current is initially so nearly 
saturated that condensation occurs immediately the motion commences and 
continues throughout the interval considered. The descending current 
produces no evaporation. Thus the parcels of air in the two currents 
follow the saturated and the dry adiabatic, respectively. Then 

o) c = r /■> %[(x-K.) * £ . 

The meaning of the various symbols above is given in the Table 

M* 

of Symbols and Abbreviations. It should be noted that M is the ratio 

of the masses of ascending and descending currents and that the cloud 

/«' M‘ 

cover due to convective activity is • As jj approaches zero the 

Slice Method reduces to the parcel method. 

Beers computational procedure consisted of combining the factor 

• M* 

R with C and making nomograms for various values of m. They were 

entered with the ratios of and J£- to find , where K is 

ft k 

Rlt, A . This multiplied by the thickness of the layer gave a measure 
of stability in arbitrary units. 

Duthie, with no further assumptions, merely by applying the 
definition of lapse rate and appropriate substitution of the hydrostatic 
equation and the equation of state, simplified equations (1), (2), and 
(3) to the following: 

I. Dry ascent-dry descent— 

(i) c. - rt&sj) (t+£)J ‘j v'&t. 

Ob- f 

0") p/£) = * j^) ■ w i . 

‘ j u~/>± j-n 
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II. Moist Ascent-Moist Descent — 

(2') C = ( r ~t- C 1 ) J ts'ot 

-=r * 

or 

( 2 ") H (£ ) 

III. Moist Ascent - M. Descent — 

(3') 0 = 4=-* [C T~.-t )*&-#£]. Snt 

or 

( 3 ") *(£)- ~+>[(*rV*(a-*Si]-'S 

Where D (£;), M J, and are dimensionless numerical 

functions which give (for a given value of ^ ) a measure of the 
stability of a layer by the Slice Method, for the conditions indi¬ 
cated, when the factor Jir is assumed constant, they vri.ll here¬ 
after be referred to merely as D, M, or W, and where no confusion 

• _ 

will result will be referred to as C. T, Tm, Td, and T* are shown 

in Figure I. The shaded portions represent equal areas. The factor 

3 

10 is inserted in order to give convenient units. 


Ay i t i-m Pt 



j ^ 


Fig. 1.- Schematic diagram of a finite layer of the atmosphere 
as represented on a pseudo-adiabatic diagram. 


S tr't 


(T+*-T.) ( l + **') . 
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It is obvious that all the variables in the above equation may¬ 
be found directly from the pseudo-adiabatic diagram and that the thick¬ 
ness of the layer is accounted for in ~f* . Consequently, computation 
is greatly simplified. Since either (1"), (2")* or (3") may be used in 
a statistical investigation, Duthie further simplified the computational 
procedure by the use of difference functions, which, when applied to the 
D function yield the IT or W function. They are 


( 4 ) = M-D = ( (l +£)-'° 3 

(5) — W — P — ( T-m- J j) - l O* 

T* 


Note that A m and w are identical except for the factor 


Now, 

(6) M&O = D + Am 

(7) W&) = J) + A„ 

A convenient computation form is shown in Figure 2. 



Fig. 2- Form for computing the stability of a layer. 
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It is obvious from the basic equations that the function D is 

always negative except when the lapse rate is super-adiabatic, and 

that the function is always positive. Their algebraic sum yields 

the value of C for the layer. In the case of moist ascent - dry 

* 

descent ^ is the same for all values of 7 ^ , whereas, for moist 

V 

ascent - moist descent it is different for each value of 7 ^ . 

The development, by Beers, of the basic equations above neglects 
frictional forces, and the term 2 cO (where F is the projection 

of the circulation path on the equatorial plane) in Bjerknes' (5) 
circulation theorem. It should be noted that absolute circulation 
and relative circulation differ only by the term 2 cO jL£ , By- 
neglecting this term we are considering the circulation in an absolute 
reference system, while the atmospheric motions observed on the earth 
are in a relative reference system. 

The following simplifying assumptions are made in all treatments 
of the Slice Method: 

1. There is no horizontal convergence or divergence. 

2. Conditions are initially barotropic. 

3. All motions are adiabatic above the surface layer. 

It is not the purpose of this paper to discuss these assumptions, which 
obviously place very severe restrictions on our atmosphere, but rather, 
to test the formulae objectively in an actual atmosphere. Petterssen 
gives a good discussion of these limitations. 

Computation by the Slice Method is begun at the Mean Convective 
Condensation Level (MCCL), which will be approximately the level of 
the bases of convective clouds, and extended layer by layer to some 
higher level. The MCCL is determined by averaging the mixing ratios 
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for the three lowest significant levels of the sounding. The point 
where this average mixing ratio line crosses the sounding on the pseudo- 
adiabatic diagram is called the MCCL. 

Beers suggests determining layers by the significant points of the 
sounding. Using Duthie’s procedure it is acceptable to take thicker 
layers when there is no significant change in lapse rate and the mois¬ 
ture distribution does not indicate the need for more layers. 

To obtain the desired measure of stability of the atmosphere Beers 

recommended summing the values of C for the various layers from the 

MCCL to the 400 millibar level. Cain, Ashley, and Champion, however, 

found that the summation to 400 millibars gave a preponderance of 

negative (stable) values. They then found that the summation to the 
o 

0 isotherm gave much better results for middle latitudes. 

The method calls for lifting the sounding 100 millibars when a 
frontal passage is expected in the forecast period, and computing the 
stability of this lifted sounding. Lifting the sounding is accomplished 
by raising each significant point of the sounding 100 millibars by 
ascending the dry adiabat on the pseudo-adiabatic diagram until satu¬ 
ration is reached then ascending along the moist adiabat. This pro¬ 
cedure is consistent with physical processes and has given good 
results. It is obvious from the orientation of dry and moist adiabats 
that a layer relatively moist at the top and dry at the bottom will 
become more stable on lifting, and vice versa. 

Other minor details of the application of the Slice Method will 
be explained in the next chapter when the need arises. 
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CHAPTER III 


COMPUTATIONAL PROCEDURES AND RESULTS 

The primary object of this investigation was to test the desirability 
of using the Slice Method of Stability analysis as an objective forecast¬ 
ing technique using single Station observations. Previous investigation 
devoted very little attention to the forecasting problem, although the 
Postgraduate School students have attempted to use the method as a fore¬ 
casting tool from time to time. In attaining this object it was hoped to 
check the conclusions of previous investigators and follow some of their 
recommendations for further investigation. 

A secondary object was to see if the value of C had a simple 
relationship with any other values which could easily be found on the 
pseudo-adiabatic diagram. 

Commencing the computation of stability at the MCCL and extending 
it to some higher level has always appeared reasonable, since we are 
primarily concerned with the energy and its changes within the cloud 
areas. The upper limit of the computation, however, has been a matter 
of considerable speculation. Cain, Ashley, and Champion found that 
using the 400 mill ibar level for middle latitudes gave a preponderance 
of negative values, but that using 0° Centigrade as an upper limit gave 
much better results. 

It was desired first of all, in this work, to determine the upper 
level of the stability computation. The 400 millibar level was not con¬ 
sidered since previous work had proved the desirability of using a lower 
level. The 500 millibar level, the 600 millibar level, and the 0° C. 
level were considered the most likely levels to give good results. 

In order to have approximately the same thickness of atmosphere for 
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consideration in every sounding it was decided to use the actual MCCL 
except in very dry air where the M3CL was above 850 millibars. In this 
case, the 850 millibar level was used. Previous experience at the Post¬ 
graduate School indicated the desirability of using this procedure. 

Using the procedure as explained in Chapter II, except as noted above, 

and using the mean temperature of each layer instead of the mean virtual 

_ • 

temperature (T*), the values of C for wet ascent - dry descent, regardless 
of moisture distribution, were computed for 7 ^ of ^ and ^ using Wash¬ 
ington, D.C. soundings. Since only a comparison of the effect of using 
different upper levels was desired, it was not necessary to lift any 
soundings for expected frontal passages. The results are shown in Table 
la. The numbers in the squares represent the number of days in each 
category. Thunderstorm verification is for the following 12-hour period. 

The tetrachloric presentation shown was found to be very simple and 
useful since a quick estimate of the accuracy of the method used can be 
made. By adding the numbers in the center and upper left squares and di¬ 
viding by the number in the upper half of the lower right square, the 
percentage of "hits" using the method is obtained. The percentage of 
"misses" can be obtained, similarly, by adding the center figure of the 
left column and the upper figure of the center column and dividing by the 
upper half of the lower right* The percentage of "hits" are shown in the 
lower half of the lower right squares. 

In this and in all subsequent data the number 0 is used as the refer- 
• * 
encp value of C. This was used because the distribution of C around 0 was 

such that moving the reference value in either direction ( * or - ) made 

practically no difference in the percentage accuracy. 

From the data in Table la it appears that the computation of C to 600 
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millibars gives greater accuracy than either the 500 millibar or the 0° C, 

/vV 3 2 

levels. The value of TvT of -jjy- gives an increase in accuracy over for 

the 600 mb* and 0° C levels, but a decrease for the 500 mb, level. In 
view of this and the fact that by using the 600 mb, and 0° 0, level more 
thunderstorms occurred with positive values of C than for 500 mb,, it 
was decided to exclude the 500 millibar level from further consideration. 
The effect of lifting these soundings for anticipated squall line 
or frontal passages is shown in Table lb. Soundings were lifted 100 

millibars, and C was computed from 900 millibars to 600 millibars using 

• ■ 

moist ascent - dry descent for all layers. This method of computing C 
for lifted soundings was used throughout this investigation and the 
results appear at least as good as those obtained by using any other 
method. 

Table lb. shows no change in the accuracy over that for the non- 
lifted computation of C to 600 mb,, but there is a slight increase over 
that for the 0° C, level. 

The previous data was based upon using moist ascent - dry descent 
for all layers regardless of moisture distribution. Previous use of the 
method indicated that this gave better results than using values of mean 
relative humidity equal to or less than 70$, as suggested by Beers, as 
the criterion for dry ascent of a layer. 
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It was desired to examine the two methods above and also the value 
of 60% mean relative humidity to determine which gave the best results. 
The 1100 soundings for July, 1947 were used. Since computation of C to 
600 mb. and an % of % had given the best results for August, it was 
decided to use those values. The results are shown in Table 2. 

Table 2 shows that all methods give very poor apparent results, 
since, had no thunderstorms been forecast for the period 1100 - 2500 for 
the entire month, the accuracy in regard to thunderstorm forecasting 
would have been 64.5%. . The best accuracy shown is 61.5%. At a glance 
there appears no choice between the 70% criterion and the exclusive use 
of wet ascent - dry descent; however, in view of the fact that not a 
single thunderstorm verified with positive values of C using the 70% 
criterion, and in view of a qualitative examination of the August 1100 
soundings it was concluded that of the two methods the 70% criterion was 
less desirable. Subsequent analysis verified this conclusion. A mean 
relative humidity value of 50% as a criterion of dry ascent, obviously, 
has little significance. 

It was next decided to compare the effect of lifting the soundings 
for expected squall line or frontal passages in the forecast period, 
with the effect of not lifting any of the soundings. The soundings were 
lifted as explained previously and C was computed for m of jz and h, 
and was summed both to the 600 mb, level and to the 0° C, level, using 
moist ascent - dry descent for all layers, regardless of moisture. 

It was found that, in most cases, the value of C increased (indi¬ 
cating greater instability), but there were some cases where lifting 

* 

decreased C. The lifted value was used, however, even when the sign of 

« 

C was reversed. The moisture distribution, of course, determines the 
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JULY 1947 


1100 Soundings 
C to 600 robs. 



/m to 




C>o 

0 

1 

1 

C<o 

11 

19 

50 


11 

20 

51 

61.5 

C7o 

2 

5 

7 

C 4.0 

9 

15 

24 


11 

20 

31. 

55 

C > o 

4 

5 

9 

C <0 

7 

15 

22 


11 

20 

51 

61.5 


Using 70% mean relative 
humidity for determining 
wet and dry layers. 


Using 50% mean relative 
humidity for determining 
wet and dry layers. 


Using wet ascent - dry 
descent for all layers 
regardless of moisture. 


Table 2. 






effect lifting will have upon the sounding* If the sounding is rela¬ 
tively dry aloft and moist in the lower layers, lifting will make it 
more unstable, 

A comparison of the effect of non-lifting and lifting for this 
period is shown in Table 3, The net result is a slight increase in 

v , /Vi' 2. 

accuracy for values of 7 * of 75 • There is no improvement for /w - /© * 
Here again* the summation of C to 600 mb, gives slightly better results 
them to the 0° isotherm. 

Table 3 shows that the highest accuracy obtained so far is obtain¬ 
ed by using the 600 mb, level and % of , and lifting the sounding 
for a squall line or frontal passage* The accuracy in respect to fore¬ 
casting thunderstorms is 76%, If no thunderstorms had been forecast 
for the 1100 - 2300 period for the entire period of July and August the 
forecaster would have been right 72*7% of the time. Obviously, the gain 
in accuracy by using the Slice method is negligible at this point. 

In view of the poor results up to this point a significance test 
was applied to the data to determine if the values of C had any signifi¬ 
cance in relation to the occurrence or non-occurrence of thunderstorms, 
or if the distribution of C and thunderstorms might have occurred by 
chance. For the data with the lowest accuracy for the two months period, 
the value of C was very definitely significant. The probability of thS 
non-occurrence of such a distribution by chance was of the order of 
*9995. 

* 

Up to this point the value of C had been, to some degree of 
accuracy, a measure of the stability of the atmosphere. Since the 
results had not yet lent themselves to objective forecasting, the problem 
became one of somehow modifying C so that it would give a dependable 
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JULY AND AUGUST 1947 


1100 Soundings 

Not lifted Lifted 



_ Bl _woT5__, 

F* 1*0 R 


C >o 

C <o 

To 600 mbs* 

c?° 

C <0 

13 

12 

25 

14 

13 

27 

/*' - Z 
A\ ~ to 

3_ 

At /O 

*1“ 

A* /O 

4 

33 

37 

3 

32 

35 

17 

45 

62 

74 

17 

45 

62 

74 

1 

L _ _ 

10 

9 

19 

n 

9 

20 

7 

36 

43 

6 

36 

42 

17 

45 

62 

74 

17 

45 

62 

76 

C >o 

CCO 

To 0° C. 

C>° 

C<o 

1 

r~___——___ 

15 

18 

33 

14 

18 

32 

2 

27 

29 

3 

27 

30 

17 

45 

62 

68 

17 

45 

62 

66 

10 

10 

20 

11 

10 

21 

n'-X 

/v\ tO 

7 

35 

42 

6 

35 

41 

17 

45 

62 

73 

17 

45 

62 

74 


Table 3. 
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value for predicting thunderstorms. 

The soundings were then examined to determine if there were some 
consistant features of the ones which had produced "misses". It was 
found that the soundings which had given positive values of C and had 
not produced thunderstorms were, in most cases, relatively dry in the 
lower layers, especially below 700 millibars. The soundings which had 
given negative values of C and were followed by thunderstorms were, in 
general, relatively moist in the lower layers. Since the MCCL was 

usually in the vicinity of the 900 mb, level, the moisture below this 

» 

level had not been considered in the computation of 0, except in 
determining the M3CL. 

A mean relative humidity for the 1000-700 mb, layer was defined 
and computed for all soundings. A measure of relative humidity, or 
"degree of saturation", seemed more desirable than a measure of absolute 
humidity. The mean relative humidity will be referred to hereafter as 
Fm for brevity. 

Fm is a weighted mean relative humidity. It is obtained by consider¬ 
ing the layers determined by the significant points of the sounding 
between 1000 and 700 mbs., and obtaining the average relative humidity 
of each layer. The average relative humidity of each layer is then 
multiplied by a factor depending upon the thickness of the layer. For 
example, if a layer is 300 millibars thick the factor is 1.0; if a 
layer is 30 mbs. thick, the factor is 0.1, etc. These products for the 
separate layers are then added to give the value of Fm. 

Using Fm as a criterion for determining the use of D, W, or M, 
stability was computed for all 1100 soundings for the following conditions: 
Case I. (1) If Fm 4 70# use D for all layers. 

(2) If Fm = 71-74# use W for all layers. 
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(3) If Fm = 75$ use M for all layers* 

(4) In the case of a lifted sounding, use the lifted value or 
that obtained from (2) or (3) above, whichever is most 
unstable. 

Case II, (1) If Fm ^ 70$ use D for all layers, 

(2) If Fm =. 71-79$ use W for all layers, 

(3) If Fm ^ 80$ use M for all layers, 

(4) In the case of a lifted sounding, use the lifted value, or 
the one obtained by (2) or (3) above, whichever is greater. 

Case III, (1) If Fm 70$ use D for all layers, 

* 

(2) In case of a lifted sounding use the lifted value of C* 

The best results for the three cases above are shown in Table 4, 
Obviously, the best results are obtained in Case III using ^ — JZ , 

Here the forecasting error in regard to thunderstorms is 11,3$ and the 
accuracy is 88,7$, This is a considerable improvement over the 76$ 
accuracy which was the best result before considering Cases I, II, and 
III, Of the 16 times that thunderstorms would have been forecast using 
this method 13 would have verified, that is, 81,2$, Of the 46 times 
that no thunderstorms would have been forecast, the forecaster would 
have been correct 42 times, that is, 91,3$ of the time. This accuracy 
is comparable with that obtained in forecasting the usual meteorological 
variables and phenomena. 

The 2300 soundings were treated separately because initial compu- 
tations of both 1100 and 2300 soundings indicated that the value of C 
for the 2300 sounding was, generally speaking, considerably more un¬ 
stable than for the 1100 sounding. It seemed to reflect the instability 
whioh preceded the sounding, and gave little indication of what would 
occur following the sounding, 
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JULY AND AUGUST 1947 


1100 Soundings 

no fZ 

1, If Fm — 70$ use D for all layers, 

2, If Fm - 71-74$ use W for all layers, 

3, If Fm =. 75$ use M for all layers, 

4, For frontal passage use lifted 
value of 5 or 3, above whichever 
is larger, 

1, If Fm = 70$ use D for all layers, 

2, If Fm = 71-79$ use W for all layers, 

3, If Fm L 80$ use M for all layers, 

4, Same as 4, above, 

1, If Fm 4 70$ use D for all layers, 

2, If Fm — 71$ use W for all layers, 

3, Use lifted sounding for passage of 
front or squall line, 

1, Same as 1, above, 

2, Same as 2, above, 

3, Same as 3, above. 

Table 4, 
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The moisture distribution of the 2300 sounding was also different 
from that of the 1100 sounding* In.general, Fm was larger, which can be 
explained partly due to cooling, but it appears that the afternoon con¬ 
vection is the primary cause of the increase in moisture, particularly 
in the upper levels* 

Only two thunderstorms occurred within 12 hours following the 2300 
sounding, consequently, no attempt was made to compare C to activity in 
the following 12 hours as had been done with the 1100 soundings. Instead 
a comparison for the following 24 hour period was made. The best results 
are shown in Table 5* 

From Table 5 it is obvious that the value of C using the 2300 sound¬ 
ing has no value as an objective forecasting tool. 

The effect upon the 2300 sounding of daytime heating was considered 
in a few cases. This was accomplished by plotting the maximum temperature - 
for the following day on the pseudo-adiabatic diagram and drawing a line 
along the dry adiabat until it intersected the sounding, C was then com¬ 
puted for this modified sounding. The effect was to increase the in¬ 
stability in most cases*. Various weight factors were applied but none 
gave any better results and will not be discussed, 

O 

It was concluded that the value of C using the 2300 sounding 
apparently has no forecasting value as presently used. If an accurate 
means of forecasting the sounding for the following day is devised, C 
may have significance in forecasting thunderstorms for that period. 

Further attempts to improve the forecasting value of the 1100 
values of C were futile. Weight factors based upon moisture and height 
were attempted, but none improved the overall results already obtained. 
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JULY AND AUGUST 1947 


2300 Soundings 
C to 600 mbs. 


K yo/T 


C>o 

13 

17 

30 

C <o 

6 

27 

33 


19 

44 

63 

63.5 


zn 

c > ° 

6 

9 

15 

c< O 

13 

35 

48 


19 

44 

63 

65 


1. Use W for all layers. 

2. Use lifted sounding for passage 
of front or squall line. 

fi ~ 3 _ 

/A lO 

1. If Fm — 70# use D for all layers* 

2, Same as 2. above. 



Table 5. 
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A time plot of C for the months of July and August is shown in 
Figures 3 ana 4. Figure 3 is a plot of both the 1100 and 2300 soundings. 

t 

0 is a summation to 600 mbs. using moist ascent - dry descent for all 
layers, regardless of moisture, and for % - ^ • Soundings were lift¬ 

ed if a squall line or frontal passage was expected in the following 12- 
hour period. The thunderstorm symbols appear at the approximate time of 
their occurrence along the time scale. 

« 

Figure 4 shows the variation with time of C as computed according to 
the method explained under Case III above, for the 1100 soundings only, 
Sinoe there is a preponderance of large negative values when using dry 
ascent - dry descent, all values of C greater than -40 were plotted 
along the -40 line. This is the data which gave the highest verification 
score of all the methods used. 

As a matter of interest it was attempted to correlate various 
quantities obtainable from the soundings as plotted on the pseudo- 
adiabatic diagram with the value of C and the occurrence or non-occur¬ 
rence of thunderstorms. The following quantities were considered: 

(1) Pressure at the 0° isotherm, (2) temperature at the 700 millibar 
level, (3) the thickness(i.e. mean temperature) of the 1000-700 mb, 
layer, (4) the pressure at the MJCL, and (5) the quantity Fm, The 
results gave no useful information and will be mentioned only briefly: 

(1) The pressure at the 0° isotherm was very spasmodic in 

O 

relation to both C and thunderstorms. Values ranged from 550 

mb, to 650 mb, with no apparent correlation. 

(2) The temperature at the 700 mb. level varied only slightly 
and had no apparent correlation with C or thunderstorms. 
Temperatures ranged from 3° to 10° C, 
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(3) The thickness of the 1000-700 mb. layer varied only 300 
feet for the entire period (between 9750 and 10050 feet) and 
for the period considered had no significance* 

(4) The pressure of the MCCL, in general, varied with the 
moisture in the lower layers of the sounding, namely, the 
quantity Fm, The relation of C and Fm is mentioned below, 

(5) Fm and C are not simply related. Although Fm is useful 

in increasing the correlation of C and thunderstorms, there is 
apparently no simple linear relationship between C and Fm. A 
scatter diagram of C and Fm is shown in Figure 5. Since no 
useful information could be obtained from a correlation coef¬ 
ficient and regression equations in view of the poor correlation, 
and in view of a limited time for this project, none were com¬ 
puted, 

A summary of the above results will be given in the next chapter. 
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CHAPTER IV 


CONCLUSIONS 

1. The Slice Method of stability analysis, while giving a reason¬ 
able measure of stability of the atmosphere, did not give a dependable 
criterion for forecasting thunderstoms using the techniques employed 
heretofore. The procedure for computing C as recommended by Beers gave 
poorer results than by using moist ascent - moist descent for all layers 
regardless of moisture. Using the latter, best results were obtained 
using -£ = i . 

2. The 2300 sounding is of no apparent value in forecasting 
thunderstorms, since in many cases it reflects the instability of the 
period preceding the sounding. 

3. Using a weighted mean relative humidity of the 1000-700 mb* 

layer as a criterion for applying moist ascent - dry descent or dry 
ascent - dry descent, and using the lifted sounding for expected squall 
line or frontal passages with £ - £ , forecasting accuracy is increased 

to an acceptable degree for short periods of time. Here the 1100 sounding 
is used for forecasting for the afternoon and evening. 

4. An accurate forecast, for a period greater than 12 hours from 

* 

the 1100 sounding, cannot be made using the value of C objectively. 

5. The period investigated is relatively short and all conclusions 
drawn from these results must necessarily be tentative until more data 
can be collected, and the methods tested on new data. It is of interest 
to note that of the three thunderstorms observed to date at Annapolis in 
1948 all would have verified using the technique mentioned in 3. above, 
and explained in Case III of Chapter III. 
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6^ There is no simple relationship between C and other variables 
found on the pseudo-adiabatic diagram. 

The following steps are presented for applying the method which 
gave the best results in this investigation, using the 1100 sounding: 
Step I, Compute a weighted mean relative humidity (Fm) of the 
1000-700 mb, layer. 

Step II, Examine the latest surface synoptic chart to determine 

if a squall line or frontal passage can be expected prior 
in the forecast period. 

Step III,If a squall line or frontal passage is expected, lift the 

sounding and compute C, using moist ascent - dry descent 

for all layers and % - jz for the layers between 900 
* 

and 600 mb. If C is negative, forecast no thunderstorms. 

♦ 

If C is positive forecast thunderstorms. 

Step IV, If no squall line or frontal passage is expected and 
Fm = 70$, forecast no thunderstorms without further 
computation. 

Step V. If Fm is greater than 70$ and no squall line or frontal 

passage is expected, compute C for moist ascent - dry 

( • 

descent from the MCCL to 600 mb. for % = . If C 

is positive forecast thunderstorms. If C is negative 
forecast no thunderstorms. 
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